Quantum-state synthesis of multi-mode bosonic fields: Preparation of arbitrary states 

of 2-D vibrational motion of trapped ions 
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We present a universal algorithm for an efficient deterministic preparation of an arbitrary two- 
mode bosonic state. In particular, we discuss in detail preparation of entangled states of a two- 
dimensional vibrational motion of a trapped ion via a sequence of laser stimulated Raman transi- 
tions. Our formalism can be generalized for multi-mode bosonic fields. We examine stability of our 
algorithm with respect to a technical noise. 
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I. INTRODUCTION 

The ability to control preparation and evolution of 
states of quantum systems opens new horizons in ex- 
perimental physics (e.g. tests of fundamental concepts 
of quantum mechanics) as well as in potential technical 
applications (such as quantum information processing). 
Recent advances in quantum optics (micromasers, cavity 
QED) Q and atomic physics (dynamics of trapped ions) 
P] have demonstrated that a macroscopic observer can 
effectively control dynamics as well as to perform a com- 
plete measurement of states of microscopic quantum sys- 
tems. In particular, preparation of an arbitrary quantum 
states of a single mode electromagnetic field in micro- 
masers and 1-dimensional vibrational motion of trapped 
ions have been discussed in literature Experimen- 
tal realizations of highly non-classical states, such as Fock 
states, squeezed states or Schrodinger-cat states of these 
single-mode bosonic fields have been reported |||| . The 
ability to synthesize an arbitrary motional state is a key 
prerequisite for a quantum measurement of an arbitrary 
motional observable of a trapped ion as proposed by Gar- 
diner, Cirac and Zoller [jio). These authors have gener- 
alized the 1-D synthesis of motional states to 2-D and 
higher dimensions. In spite of the conceptual elegance, 
the method proposed by Gardiner, Cirac and Zoller is dif- 
ficult to implement in general. The problem is that the 
number of laser operations required for a preparation of 
a two-mode target state 



target , 



of the form 



I* 



M ma x 



N max 



target) 



m— n—0 



Q mn \m, n) 



(1.1) 



depends exponentially on the dimensionality of the sub- 
space of the Fock space in which the target state is em- 
bedded. In particular, if N max — M max then the num- 
ber of necessary operations is proportional to 2M max x 
2 m ™ . This exponential dependence restricts applica- 
bility of the proposed method (see below). A novel ap- 
proach which overcomes this obstacle was introduces very 
recently by Kneer and Law jnj. The authors have con- 
sidered a photon-number-dependent interaction which is 



induced by a detuned standing-wave laser field in y di- 
rection. It keeps on resonance only transitions between 
Fock states \m, n) with a fixed number n — n y of trap 
quanta in y direction while applying sequentially the 1-D 
preparation scheme |8| in the x direction for each partic- 
ular n = n y . 

In this paper we propose an alternative universal algo- 
rithm for a construction (synthesis) of an arbitrary quan- 
tum state of the two-mode bosonic field which can be 
straightforwardly generalized for the preparation of an 
arbitrary multi-mode system. In our algorithm the num- 
ber of operations required for a preparation of the two- 
mode state ([0]) grows only polynomially as a function of 
Nmax and M max . In particular, if N max = M max , then 
the number of operations is proportional to 8M „ 



I 2 . 

max 



To make our discussion as close as possible to exper- 
imental realization we consider the preparation of the 
two-dimensional vibrational motion of a trapped ion p2[ . 
The choice of the system is motivated by the fact that 
dissipative effects can be significantly suppressed in ion 
traps which is important for a deterministic engineering 
of quantum states. Moreover, the quantized vibration 
motion of a trapped ion can be effectively controlled by a 
proper sequence of laser pulses tuned either to the atomic 
electronic transition or to resolved vibrational sidebands 



The simplest quantum state preparation is represented 
by a simple unitary evolution of an input state of the sys- 
tem governed by a specific (generally, nonlinear) Hamil- 
tonian. Obviously, the fixed Hamiltonian restricts the 
family of target states which can be "generated" from 
available inputs. Another way how to prepare states of 
the system of interest (e.g, single-mode bosonic field) is 
to consider quantum interaction between this system and 
an adjoint quantum system (e.g. fermions). Dynamics of 
these two systems is governed by a specific interaction 
Hamiltonian. The desired quantum state engineering is 
then achieved by conditional measurements performed on 
the adjoint quantum system Q. This non-unitary selec- 
tion of specific quantum trajectories allows to synthesize 



1 



essentially all quantum states of the system under con- 
sideration. But there is prize to pay - the probability of 
the outcome of the given conditional process can be ex- 
tremely small. To overcome this problem one may con- 
sider a sequence of interactions between the original and 
the adjoint systems. These interactions (channels) are 
governed by different Hamiltonians with just one chan- 
nel turned on at a given time. Coupling constants in 
these Hamiltonians and times of duration of given in- 
teractions are in this case free parameters which can be 
appropriately tuned so that at the output the system of 
interest is disentangled from the adjoint system (see be- 
low) and is prepared in the desired state. This approach 
has been recently utilized by Law and Eberly || who 
have shown that by a suitable switching between two 
channels of the atom-field interaction one can generate 
an arbitrary state of the single-mode cavity field. It is 
important to note that quantum states of 1-D vibrational 
motion of trapped ions are experimentally created also 
via a sequence of laser pulses tuned either to an elec- 
tronic transition or an appropriate vibrational sideband 
that is, sequential switching of different interaction 
channels is used. Generalizations of the 1-D vibrational 
quantum state synthesis to 2-D and more dimensions 
have been discussed recently [ ^0|JT]] ]. In what follows we 
introduce an universal scheme which enables determin- 
istic preparation of an arbitrary two-mode target state 
via an appropriate switching between laser stimulated 
Raman transitions described by single and two-mode in- 
teraction Hamiltonians. The method can be generalized, 
e.g., for a trapped ion in a 3-D trap potential and other 
multi-mode bosonic systems. Even for higher dimensions 
the number of preparation operations scales polynomially 
with the dimensionality of the subspace of the Fock space 
in which the target state is embedded. 

The paper is organized as follows: in Section II we 
introduce our algorithm for 2-D quantum-state synthe- 
sis. In Section III we discuss a physical realization of the 
preparation scheme. In Section IV we examine the stabil- 
ity of our algorithm with respect to a technical noise. Fi- 
delity between outputs and specific target states is eval- 
uated. We finish our paper with conclusions. 



II. QUANTUM-STATE SYNTHESIS 

Without any loss of generality let us assume that we 
want to generate an arbitrary two-mode sta te g iven as a 
finite superposition of number (Fock) states (1.1). To de- 
scribe the preparation algorithm we firstly split the whole 
Hilbert space into an appropriate set of finite-dimensional 



subspaces labeled by a specific quantum number. Then 
we introduce two sorts of interaction channels. The 
first set of Hamiltonians will be responsible for dynam- 
ics within a given subspace; while the second one will 
realize a "transfer" of probability amplitudes between 
different subspaces. In particular, let us consider vibra- 
tional states of a trapped ion, confined in 2-D harmonic 
potential. Excitations of two vibrational modes are de- 
scribed by creation and annihilation operators d M and aj, 
(p = x,y). The Hilbert space H V ib of all two-mode vibra- 
tional states can be divided to subspaces with constant 
total number of vibrational quanta, i.e., 7i v a = 
J = 0, 1, ... , oo, where Tlj is spanned by two-mode Fock 
states | J, 0), | J — 1, 1), . . . , |0, J). The subspaces Hj are 
dynamically independent (invariant) for those Hamiltoni- 
ans which do not change the total number of vibrational 
quanta J — al,a x + a y a y . However, to realize controlled 
manipulations with states within a given subspace TLj we 
need an interaction of the 2-D bosonic field with another 
quantum system. In the case of trapped ions it is natural 
to assume the coupling to internal energy levels of ions. 
In particular, we use three internal electronic states \i) 
(in A or 3 configurations) with energies fiUi (i = a, 6, c) 
which form the basis of the Hilbert subspace Tti n . Due 
to quantum interaction between the vibrational and in- 
ternal degrees of freedom the total state vector of the 
composed system with the given maximum number of 
vibrational quanta J m ax = M max + N max at time t reads 

M max N max _^ 

!*(*))= J2 Qm,n;i(t)\m,n)®\i) 

m—0 ii—O i—a 7 b,c 

= EE E Qk,J-k;i(t)\k,J-k)®\i), (2.1) 

which reflects quantum-mechanical entanglement be- 
tween the two subsystems. We assume that initially the 
ion is in an internal state |a) and the vibrational motion is 
cooled to the ground state |0, 0). That is, the initial state 
vector of the composed system \*S?(t = 0)) = |0,0) <g> \a) 
is factorized. Our task is to find an unitary evolution 
U\ such that at some time t = T the state vector 
( |2.lD can again be factorized, while the vibrational state 
of the ion is described by the target vector ([__]), i.e. 

|*(t = T)) = #t| 0) 0) ® \a) = ^target) f |o>. 

In what follows we prove that the unitary transforma- 
tion U can be represented by a sequence of five "elemen- 
tary" unitary transformations (p = 1,...,5) which 
act in the product Hilbert space H. V ib <S> Ti-in and which 
correspond to five interaction channels associated with 
interaction Hamiltonians H^: 



U = A { J] CjBj^A.,}- Aj = U\y m ult 1Xa 



M 3) 



(i) 



|2,J-2,6)^|l,J-l,a>^|l,J-l,6)^|0,J,o>' 



B, 



j=i 
= U, 



(2) 



M 4) 



(2) 



| J,0,c> |J-1,1,6> |J-l,l,c) 



.J7, (4) 



M 2) 



u. 



(4) 



(2) 



|l,J-l,b) w |l,J-l,c)^|0,Xb)HO,J,c)' 



c, = u, 



(5) 

I J,0,b) 



(2.2) 
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FIG. 1. Recursive "de-evolution" algorithm: (a) Vectors |m,n) ® \i) are represented as points of a lattice m,n with internal 
levels i = a,b,c shown in "ovals" . Components which enter (at this stage of "de-evolution" ) the state vector ( |2.l[ ) are indicated 
by •. Dashed lines connect basis vectors from subspaces with constant number of trap quanta J. (b) The action of the operator 

a) to the 



Aj for J 



>. The "elementary" transformation j ^ transfers completely the population of the state |0, J) 



state |1, J — 1) ® |6) as indicated by solid arrow. The atomic transition scheme of the corresponding stimulated Raman process 
p = 3 is shown in the inset. Next 1 ^ transfers the population from |1, J — 1) <g> |6) to |1, J — 1) ® \a) by means of the 

interaction channel p = 1. The dashed arrows indicate simultaneous transitions which lead to a change of the state vector ( |2.l[ ). 
The application of the sequence Aj "shrinks" the whole population from the given subspace Hj Cg> Tii n to the component state 
[1, J — 1) ® a), (c) Specific preparation of the subspace Hj_ 1 (8>Tlin by the operator Bj_ 1 - The operator l ^ transfers 

transfers the population from |0, J — 1) Cg> |i>) to 



(4) 



the population of the state |0, j — 1) ® |c) to |0, J — 1) CS> |6), while L/ Q - 



J-1,6) 



[1, J — 2) (g> |6). The inset shows the schemes of the utilized stimulated Raman transitions p = 2,4. Finally, the sequence Bj_ 1 
leaves on ®TLi n only the component state | J— 1,0} £g> |6) and those states with the internal level |a) contributing to (2.1). 

(d) The operator Cj = Uff n n ^ transfers the population of the component state J, 0) ® \a) to the state |J — 1,0) ® |i>). The 

a) into |0,0) ® \a). 



J J "\J,0,a) 

procedures (b)-(c) are recursively repeated to "de-evolve" the initial state \iptarget) 



The operators Aj, Bj, Cj are build up from "ele- 
mentary" unitary transformations ufj?j_ k * where su- 
perscript p determines the interaction channel and sub- 
script indicates that parameters of the interaction are 
adjusted to transfer completely the current population of 
the component state |fc, J — k) eg) \i) to the "neighbor- 
ing" state according to the given type of interaction (see 
below). 

The set of the "elementary" unitary transformations 
is sufficient (even though not unique) to evolve the 
target state (1.1) into the vacuum and vice versa. To 
specify them we will analyze the "de-evolution" of the 



target state into the vacuum state. We will perform this 
de-evolution via a systematic "transfer" of population 
from state vectors with higher to smaller number of vi- 
brational quanta. That is, we choose the operators 
so that at a given step of the unitary de-evolution a par- 
ticular probability amplitude Qk,j-k-i is made equal to 
zero. In Fig.l. we visualize the "direction" of the action 



of operators U, 



(p) 

\k,J-k,i) ■ 



Namely, in Fig. la we represent 



state vectors \m,n) (g> \i) as "ovals" located at the points 
of the lattice marked by parameters m, n associated with 
vibrational state of an ion, while its internal state repre- 
sented by a simple energy level diagram [population of a 
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given internal level is indicated by a bullet (•)]. Dashed 
lines indicate the subspaces Hj ® Hi n labeled by a con- 
stant number of trap quanta J. 

Let us s tart the "de-evolution" procedure of the target 
state (1.1) into vacuum. After time t = T — t the vector 
describing the vibra tional and internal state of the ion 
is given by Eq.(2.1). We assume that the de-evolution 
is performed in such way that at this moment the state 
vector (2.1) is composed only of state vectors from sub- 
spaces Hj ® Hin with J < J (< Jmax)- Now we apply 
three sequences of operations described the operators Aj, 
Bj_ 1 , and Cj, respectively [see Eq.(|2.2|)]. 

The action of the operator Aj is illustrated in Fig. lb. 
It describes a process in which the first "elementary" 
transformation £/? - ^ transfers completely the popula- 
tion of the state |0, J)®\a) to the state |1, J— 1)(£>|&) (see 
the solid arrow in Fig. lb; dashed arrows indicate simul- 
taneous transitions which are not controlled at the given 
stage). This process can be realized by irradiating the ion 
with two external laser fields with tunable frequencies io x 
and uj y in x and y directions. Adjusting the resonance 



conditions cu v 



LUb — uj a + v x — v y (= lu 3 ) the stimu- 



lated Raman transition [see the inset in Fig. lb, p — 3] can 
be described in the Lamb-Dicke regime by the effective 
interaction Hamiltonian (we discuss physical conditions 
under which this and all subsequent Hamiltonians can be 
justified in Section III): 



H 3 =g 3 aiay\b){a\e- i ^ t 



g 3 a x a y \a 



t l-)(6|e^ 3t ; 



(2.3) 



with the corresponding time evolution operator = 
exp(— iH^H) setting h — 1. If the interaction constant 
.93 = | <73 1 e l6>3 and the duration of the interaction t are 
chosen to fulfill the condition 



^ 3 Qk,j-k-a cos(| ff3 |*V(fc + ±)(J- k)) + 
)k+i,J-k-i;b sm(\g 3 \ty/(k + 1)( J - k)) = 0, 



(2.4) 



then the population of the component state \k, J—k)®\a) 
is completely transferred to the state |fc+l, J— k— l)®|i>). 
(for instance, in the situation described by Fig. lb the 
probability amplitude Q j. a becomes equal to zero). 

Once this is done, then the transformation U,}\ ,. is 

|1,J— 1,6) 

turned on. In this process the population of the com- 
ponent state |1, J — 1) ® \b) is completely transferred to 
the state |1, J— 1) (8 \a). The corresponding interaction 
channel is described by the Hamiltonian 



ffi=fli|6)(o|e- 



■g\\b) (a\e 



iuj\t 



(2.5) 



Here we assume lasers to be tuned to the electronic tran- 
sition, i.e., uj y — lu x — LUf, — io a (= u>i ). To cancel the 
term Qk,j-k-b m the state vector (2.1), the interaction 
constants (|<7i|t, 0i) have to be chosen to satisfy the con- 
dition 

Qk,J-k;a sin(lfil*) + i e ~ l9l Qk,j-k;b cos(lsnjt) = 0. (2.6) 



In particular, in Fig. lb the operator 



(i) 

|l,J-l,o) 



cancels 



J-l:b 



for J = 3. The successive action of the "ele- 
mentary" transformat ion s and which form the 
operator Aj [see Eq.(2.2)] finally "shrinks" the popula- 
tion of the subspace Hj®Hi n to a single state | J, 0)£g>|a). 

At this stage we start the process of cancellation of 
the contribution of component states \k, J — 1 — fc) (8) |6) 
and \k, J — 1 — k) ® \c) in the "neighboring" subspace 
Hj_ j® Hin (see Fig.lc). This intermediate procedure is 
required to prevent a reverse transfer of population from 



H 



Mi 



to Hji 



)TLin (see below and Fig. Id). For this 



purpose the operator Bj_ 1 is constructed from "elemen- 



tary" operations and [see Eq.(2.2)]. Namely, 

" (2) 

the operator U, = , describes the transfer of popula- 

|0,J— l.c) 

tion of the state |0, J - 1) <g> |c) to |0, J ~ 1) (g) |6). This 
transfer can be achieved with the help of lasers pulses 
tuned to the electronic transition between the levels \b) 
and |c). The corresponding interaction Hamiltonian 
and t he r esonance condition are analogous to those for 
.ffW ( |2.5| ) [we have to replace only b — > c and a — > b]. 
Further the operation U^ 4 \ , , cancels a contribution of 

|0,J— l,o) 



the state |0, J — 1) ® \b) to the state vector (2.1). This 



interaction channel is described by the Hamiltonian H 1 - 4 ^ 
which is obtained from (^^) by the substitution b — > c 
and a — > b. We see that the operator Bj_ 1 acts like 

Aj but instead of the stimulated Raman processes be- 
tween a <-> b the transitions b *-» c are utilized. As the 
result of the action of the operator Bj_ 1 in the subspace 



n 



.7-1 



Hi 



only the component states with the inter- 
nal level | a) and the state | J — 1,0) ® \b) have non-zero 
amplitudes. 

After the action of the operators Aj and Bj_ x the uni- 
tary operator Cj is utilized to transfer the population of 
the component state | J, 0)<£>|a) to the state | J— 1, 0)® |6) 

(see Fig. Id for J — 3). The transformation , which 

\Jfi,a) 

performs transitions between subspaces with the number 
of trap quanta differed by one is realized by the process 
described by a single mode interaction Hamiltonian 



H 5 = grMb)^- 1 ^ + gtal\a){b\i 



(2.7) 



with the resonance condition uj y — u) x = Ub — uj a — v x 
(= W5). The parameters of this interaction channel are 
determined by the constraint 

ie ifls Qj,o ; aCos(|s 5 |*v7) +Qj-i,0;bsin(\g 5 \tVJ) = 0, (2.8; 
for J = J. 

As the result of the action of the operators Aj, Bj,, 
and Cj, all coefficients Q k j_ k . i in (2.1) are equal to 



zero. Moreover, the situation before (Fig. la) and af- 
ter (Fig. Id) the action of these operators is the same, 
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except we just "moved" from the subspace Ttj ® TUn to 
H 



j-i 



Ttin with the number of trap quanta decreased 
by one. This means that the procedure can be recur- 
sively repeated. The given solution of the "de-evolution" 
gives immediately the recipe for the creation of the target 
vibrational state from the two-mode vacuum (with elec- 
tronic level | a)): We have to change by tt the phase shift 
between external laser fields and repeat the sequence in 
the opposite order, i.e. we apply W on |0, 0) ® \a). In 
this way we can synthesize an arbitrary state. The num- 
ber of "elementary" operations involved in this process 
is proportional to 2J£ lax , which is important for an ex- 
perimental realization of the proposed scheme, i.e., the 
number of necessary operations increases only polynomi- 
ally with the increase of the size of the Hilbcrt space in 
which the target state is embedded. 

In the following section a physical implementation of 
the preparation scheme is discussed in more details. 



III. REALIZATION OF INTERACTION 
CHANNELS 

The Hamiltonians which are eligible for the synthe- 
sis of two-mode bosonic states has been discussed re- 
cently by several authors Some of the propos- 
als are based on laser stimulated dipole transitions |Tc| ] 
and phonon-number-dependent interaction via detuned 
standing-wave JTlj] . For our purposes we have utilized 
the stimulated Raman processes discusses in detail by 
Steinbach, Twamley and Knight (Rj. In this section fol- 
lowing Ref. |]l3f we briefly derive the Hamiltonians which 
are used in our algorithm [see Eqs.(|2.3|),(p^),(2.7)] and 
discuss the range of their applicability. Let us consider a 
trapped ion confined in a 2D harmonic potential charac- 
terized by the trap frequencies v x and v y in two orthog- 
onal directions x and y. The ion is irradiated along the 
x and y axes by two external laser fields with frequencies 
uj x , ojy and wave vectors k Xl k y . The laser fields stimulate 



Raman transitions between two internal energy levels \a) 
and | b) via an auxiliary electronic level |r) which is far off 
resonance. For concreteness, we consider A-configuration 
with the upper level |r) as outlined in the inset of Fig. lb. 
The interaction Hamiltonian for the system under con- 
sideration can be written in the dipole and rotating-wave 
approximation (RWA at laser frequencies) in the form: 

H in t =g* x S k **-»^\r){b\ +g x e-^~^\b){r\ 

+ g* y e l{k ^- bJ y )t \r){a\ + gye-^-^la}^]. (3.1) 

The coupling constant g x (g y ) is proportional to the in- 
tensity of the laser field in x (y) direction and the dipole 
moment of the electronic transition \b) <-> |?') (\a) «-> \r)). 
The upper off-resonant level |r) can be adiabatically 
eliminated provided that A x ,A y ^ g a ,<?b,|A x — A y |, 
where laser detunings for dipole transitions \b) «-> |r) 
and \a) «-> |r) are denoted as A^. = (uj r — cl>{,) — ui x and 
A y = (uj r — u) a ) — Wy, respectively. After adiabatic elimi- 
nation the effective interaction Hamiltonians for the stim- 
ulated Raman transition \a) <-> \b) reads jl3| 

Ht{ f) = fe-^^M-ie^b^e^b)^ 

+ g e i ( w »- w -J*D ar (ie B )D v (-ie y )|o)<6|. (3.2) 



Here D q (ie q ) = c 



ie q (a q +a q ) 



^ik q q 



is the displacement 



operator (q = x,y); the Lamb-Dicke parameter e q is de- 
fined as e q = h 2 k q I '{2mhv q ) and the effective interaction 
constant g — g x g y (-^-+-^-)- Further, we assume that the 
energies of the electronic levels \a) and \b) are redefined 
to include Stark shifts due to the adiabatic elimination 
of the off-resonant energy level |r) 

In the interaction picture the effective interaction 



~ (eff) 



in I 



e iH tft{eff) e -iH t can be ex _ 



Hamiltonian H 

pressed as [the free Hamiltonian Hq induces transforma- 



tions a„ 



a a e 



|a)(6| -> laXble-*^*-"")*] 



■x ^ff) 



9 e 



{-ie x ) 



k+m (ie y ) l+n 



-it[A.+(k—m)v x + (n—l)v y ] *t m 



k\l\m\n\ 



al m a k x a y l a n y \b)(a\ 



(3.3) 



The resonant terms in the expansion (3.3) which con- 
tribute dominantly to the resulting effective Hamiltonian 
can be selected by an appropriate choice of laser frequen- 
cies. If off-resonant processes are oscillating with suffi- 
ciently high frequencies they can be eliminated applying 
the second RWA at trap frequencies. 

In particular, tuning lasers to the first red sidebands 
A = A x — A y = v x — v y with incommensurate trap fre- 
quencies only the resonant terms with k_— m = 1 and 
I - 



1 are retained in the expansion (3.3) 



- (3) 



H int = g*e x e v alT{ala x , a\a v )a v \b) {a\ + h.c. (3.4) 



where 



T 



(-1) 



k+l 2k 21 



kj 



(fc + l)!fc!(Z + l)!Z 



-tf k a k a) l a l 

I u x u x u y u y 



(3.5) 



In the Lamb-Dicke regime e x ,e y <C 1 the operator T 
is close to the unity operator and the Hamiltonian (3.4) 
which is written in the interaction picture acquires in the 
Schrodinger pict ure exactly the form of the two-mode 
Hamiltonian H3 (2.3). With proper laser tunings we can 
design within the Lamb-Dicke limit also remaining in- 
teraction Hamiltonians required for the quantum state 
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synthesis described in S ecti on II. In particular, retaining 
only resonant terms in (3.3) for the laser-stimulated Ra- 
man process with A = the interaction Hamiltonian Hi 
(2.5) is obtained in the Schrodinger picture. The process 
with A = — v x is described by the one-mode interaction 
Hamiltonian H$ (2.7). Let us also notice that beyond the 
Lamb-Dicke limit this process is analogous to a nonlin- 
ear Jaynes-Cummings dynamics discussed by Vogel and 
de Matos Filho @. 

Considered approximations impose limitations on the 
applicability of the inte raction Hamiltonians. Driven 
electronic transition 



( |2.5| ) and one-mode interaction (2.7) 
have been considered already for the 1-D quantum state 
synthesis A new tool in our approach represents 

the two-mode interaction Hamiltonian H3 (|2.3|). Com- 
prehensive analysis of the limitations for ( |2 .3|) was done 
in Ref. Jl3fl . It turns out that the most subtle poi nt i s 
the second RWA at trap frequencies [see Eqs.([L3|),(3.4)1 
which imposes restrictions fl3|] 



\g\e x e v max(N max , M rnax ) < mhx(u x , u y ), 
min(^ x , u y ) 



(3.6) 



max(z/ x , u y ) 



> 5. 



The trap anisotropy, i.e., incommensurate trap frequen- 
cies, is required to avoid additional resonances in ( |3.3| ). 
Numerical simulations in Ref. |l3| demonstrated that it is 
experimentally feasible to operate the considered Hamil- 
tonians within the Lamb-Dicke limit. 

It is worth noticing that already the nonlinear form 
(3.4) corresponding to the two-mode interaction Hamil- 
tonian ( p. 3[ ) outside of the Lamb-Dicke regime al- 
lows us to adopt the proposed algorithm for quan- 
tum state synthesis. The main difference consists in 
the form of generalized Rabi frequencies. In partic- 
ular, the matrix element of the interaction Hamilto- 
nian (3.4) in the Lamb-Dicke regime e q -C 1 [i 



1. 



i.e. 

m 4 



(2.3) in the interaction picture] reads (b,n 

~ (3) 

l\H int \m, n, a) = g^y (m -I- l)n while beyond the Lamb- 
Dicke regime the corresponding matrix element is given 

as g 3 e- J1 ^ JL Ll n (e 2 x )Ll l _ 1 {e 2 y )/y / (m+ Tjn where L\ n is 
the associated Laguerre polynomia l. In s erti n g th e "non- 
linear" Rabi frequencies into Eqs.(2.4),( |2.6| ),(2.5) which 
determine the choice of the interaction constants and 
switching ti mes, the same sequence of elementary trans- 
formations (2.2) could be applied for preparation of a 
given state even outside of the Lamb-Dicke regime. On 
the other hand, the conditions for the applicability of the 



second RWA at trap frequencies to obtain (3.4) from (3.3) 



are not so transparent as res trictions for the RWA within 
the Lamb-Dicke regime (3.6). This problem goes beyond 
the scope of the present paper in which we demonstrates 
the potential of our algorithm in the Lamb-Dicke regime. 



IV. STABILITY OF SYNTHESIS 

Our algorithm works ideally when there are essen- 
tially no dissipations in the system. It is the reason 
why we have consider effectively dissipation-free dynam- 
ics of trapped ions. Nevertheless, we have to stress that 
the ideal synthesis of quantum motional states assumes 
also perfect control of interaction constants and switch- 
ing times of particular interaction channels. Namely, 
the values of g p t foun d via th e "de-evolution" procedure 
[i.e., solutions of Eqs.(2.4,2.6,[2.8])] have to be controlled 
as precisely as possible. In practice one cannot avoid 
some level of "technical" noise, for example, due to im- 
perfect timing of switching between interaction channels. 
Therefore in what follows we will study the stability of 
the presented algorithm with respect to the "technical" 
noise. This noise is simulated as random fluctuations 
of the ideal values g p t. In particular, fluctuations are 
equally distributed aro und the ide al (complex) values g p t 
[i.e., solutions of Eqs.(2.4,2.6, 2.8)] within a fixed inter- 
val (1 + i)5. In general, the output state vectors \^s) 
of the composed system which are prepared in the pres- 
ence of "technical" fluctuations take the form (|2.l|), i.e., 
the vibrational and internal degrees of freedom are not 
disentangled. The Fig. 2 shows the fidelity / of the "im- 
perfect" output states \^s) with respect to the desired 



state vector \^ target) ® |a). The fidelity is defined as the 
averaged squared scalar product of particular realizations 
with \^ target) ® \a), i.e., 



f = ms\*target)\a)\ 2 ))s. 



(4.1) 



In our simulations we have performed averaging over 
100 runs of state-synthesis sequences. In these runs 
each value q p t associated with a given elementary op- 
eration acquires a random fluctuation within the 
interval (l + i)8. We have considered two different target 
states, two-mode cat-like state \S& C at) — N ca t(\ct)\a) + 
I — a) I — a)) and two-mode correlated state \^ C orr) — 
e -M /2 ^=| TOj m ), with a = 2. For this value of the 
amplitude these two states have approximately the same 
mean number of vibrational quanta n — n x 4 n y ~ 8.0. 
We have considered two cases when M max — 12 and 
M max = 20 (here N max = M max ). 

Fig. 2 clearly indicates the fact, that larger the value 
of M max more pronounced is the role of fluctuations 
(compare solid and dashed lines which correspond to 
M max = 12 and M max = 20, respectively). This ob- 
servation is easy to explain: The total number of oper- 
ations in our algorithm is proportional to 8M^ ax (for 
M max = N max ) which means that the case M max = 20 
requires almost three times more operations compared to 



the case with M„ 



12. The noise is accumulated as a 



function of elementary operations, therefore to improve 
the fidelity of the preparation process it is important to 
choose carefully M max . To be specific, for a given e we 
have to choose the minimal values of M max and N max 
such that 
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N„ 



EEiQ™i 2 - E E \Q m n\ 2 <s. 



(4.2) 



m—0 n—Q 



m— n— 



We can also use Fig. 2 to illustrate the fact that for a 
given value of M max the fidelity of the preparation may 
depend on the target state. To be specific, we have found 
that smaller the number of nonzero amplitudes Q mn in 



the target state (1.1) the higher the fidelity is for a given 



value of the range of fluctuations S. This behavior can be 
rather surprising as for synthesis of two states "localized" 
within the same region of the vibrational "lattice" m, n 
we need comparable number of elementary operations. 
Nevertheless, each nonzero Q mn after the synthesis is bi- 
ased by more or less the same error and, consequently, 
the states with smaller number of nonzero amplitudes 
Qmn are less sensitive to fluctuations. As already stated, 
this observation concerns only states "localized" within 
the same region of the vibrational lattice m, n and with 
a comparable mean number of vibrational quanta. 




states with relatively small number of vibrational quanta. 
This means that not only the number of operations is ex- 
ponentially growing but also the number of vibrational 
quanta during the preparation procedure may transiently 
exponentially increase. 

We stress that our algorithm is associated with manip- 
ulations only within the original subspace of the Hilbcrt 
space specified by M max and N max . Moreover, the num- 
ber of operations grows only polynomially as %M"^ lax . 
The dimension of the two-mode Fock subspace from 
which the component states are does not increase dur- 
ing the preparation procedure. This great reduction of 
number of operations is due to the fact that we have 
employed the third atomic level \c) in the preparation 
procedure. Very recently a new 2-D preparation scheme 
was introduced by Kneer and Law pT| . The standing - 
wave laser field in y direction induces a photon-number 
dependent interaction enabling thus to use 1-D schemes 
selectively for particular subspaces with constant number 
of trap quanta in y direction. The number of operations 
scales also polynomially as 2M% lax with only two elec- 
tronic levels involved. 



0.9 



u 0.8 



0.7 



0.6 




0.0 0.005 



0.01 

6 



0.015 0.02 



FIG. 2. We plot fidelity given by Eq.(^) as a function 
of the range of fluctuations 5. We consider two different tar- 
get states, two-mode cat-like state |$ C at) (denoted by A) 
and two-mode correlated state \^ CO rr) (denoted by rj) with 
a = 2. We have considered two cases when M max = 12 (solid 
lines) and M ma x = 20 (dashed lines). Comparing two solid 
lines we see that smaller i s the number of nonzero amplitudes 
Qmn in the target state (1.1) the larger is the fidelity for a 
given value of 8. 



V. CONCLUSION 

In this paper we have presented a universal algorithm 
for an efficient deterministic preparation of an arbitrary 
two-mode bosonic state. We have adopted this algorithm 
as a computer program The proposed method can 
be generalized to 3-D trapping potential and three-mode 
vibrational states. In this case one would need four in- 
ternal electronic levels and nine interaction channels cou- 
pled to the 3-D vibrational field. It can be shown that the 
number of operations required for quantum state synthe- 
sis scales polynomially (~ M^ ia:r ). Further generalization 
to multi-mode fields is possible fl5fl. 
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Finally, we briefly compare our algorithm with the one 
proposed by Gardiner, Cirac and Zoller jl0| in which 
the number of operations in the preparation sequence 
is growing exponentially as 2M max x 2 max . Therefore 
fluctuations for M max large enough cause insurmount- 
able problem. There is also another problem with this 
procedure. Namely, Gardiner, Cirac, and Zoller have 
utilized only two internal atomic levels (\a), \b)) which 
results in the fact that their algorithm is based on ma- 
nipulations with vibrational states which are out of the 
original Hilbert space specified by the cutoffs M max and 
N m ax- In other words, manipulations with highly ex- 
cited vibrational states are required for construction of 
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